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Abstract 

The super diffusion behavior, i.e. < x 2 (t) >~ t 2u , with v > 1/2, in general is not 
completely characherized by a unique exponent. We study some systems exhibiting 
strong anomalous diffusion, i.e. < >~ t qu ^ where u(2) > 1/2 and qv(q) is 

not a linear function of q. This feature is different from the weak superdiffusion 
regime, i.e. v(q) = const > 1/2, as in random shear flows. 

The strong anomalous diffusion can be generated by nontrivial chaotic dynamics, 
e.g. Lagrangian motion in 2d time-dependent incompressible velocity fields, 2d 
symplectic maps and Id intermittent maps. Typically the function qv{q) is piecewise 
linear. This corresponds to two mechanisms: a weak anomalous diffusion for the 
typical events and a ballistic transport for the rare excursions. 

In order to have strong anomalous diffusion one needs a violation of the hypoth- 
esis of the central limit theorem, this happens only in a very narrow region of the 
control parameters space. 

In the presence of the strong anomalous diffusion one does not have a unique 
exponent and therefore one has the failure of the usual scaling of the probability 
distribution, i.e. P(x,t) = t~ u F(x/t"). This implies that the effective equation at 
large scale and long time for P(x,t), cannot obey neither the usual Fick equation 
nor other linear equations involving temporal and/or spatial fractional derivatives. 

PACS number(s): 05.45.+b, 05.60. +w: 
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1 Introduction 



The transport of mass or heat from particles passively advected by a given velocity field is 
a problem of considerable practical and theoretical interest. The coupling between advec- 
tion and molecular diffusitivy results in a great variety of different fields, as geophysics, 
chemical engineering and disordered media |2j . One can have highly nontrivial behavior 
even in the presence of simple laminar velocity fields, e.g. very large diffusion coefficients 
due to the combined effects of the molecular diffusitivy and some features of the velocity 
fields 11. 

Taking into account the molecular diffusion, the Lagrangian motion of a test particle 
is described by the Langevin equation : 

x = u(x,t) + $ > (1) 

where u(x,t) is the Eulerian velocity field at the position x and the time t, and £ is a 
Gaussian white noise with zero mean and correlation function 

mW)) = ^5, l3 5{t-t') ) (2) 

the coefficient D being the molecular diffusivity. If G is the density of tracers, the 
Fokker-Planck equation associated to (|l]) is 

d t Q + d(u G) = D Q d 2 e. (3) 

For times much larger than the typical time of u, the large-scale density field (9) (i.e. the 
field 9 averaged over a volume of linear dimension much larger than the typical length of 
the velocity field u) obeys a standard diffusion Fick equation: 

dt{Q) = Dpl X] {Q) i.j 1 d. (4) 

All the (often nontrivial) effects due to the presence of the velocity field are in the eddy 
diffusion coefficient Df-. Of course if the Eq. ([|) holds then one has ((x(t) — x(0)) 2 ) ~ 
2 Df x t. In practice at large time the test particle behaves as a Brownian particle. 

The above scenario is the typical one. Nevertheless there exist cases where anomalous 
diffusion is observed, i.e. ((x(t) — x(0)) 2 ) ~ t 2u with v ^ 1/2. The case when v < 1/2 
corresponds to subdiffusion while for v > 1/2 one has super diffusion. Trapping is the 
basic mechanism leading to subdiffusion and it occurs only in compressible fields. A 
remarkable example of subdiffusion, ({x(t) — x(0)) 2 ) ~ [ln(t)] 4 is the random walk in Id 
random potential 0]. 

If the velocity field is incompressible and the molecular diffusivity is non zero, either 
standard diffusion or super diffusion takes place. 
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In the last two decades many authors have studied systems with anomalous (non 
Gaussian) diffusion. It is practically impossible to cite the huge literature. We just men- 
tion, among the many, the contribution of Geisel and coworkers || for the Id intermittent 
maps, Ishizaki et al |J for the symplectic maps, Zaslavky et al j7| for the lagragian motion 
in q-flows [§, ||, Majda & Avellaneda |T(J for the random shear flows, the experimental 
study of Solomon at al. fill on the superdiffusion in an anular tank, the symbolic dynam- 



ics approach of Misguich et al. |nj for the subdiffusive behavior in a stochastic layer and 
the anomalous diffusion of magnetic lines in the 3d turbulence by Zimbardo et al. [13| . 



The aim of this paper is a detailed study of the anomalous diffusion in simple veloc- 
ity fields and discrete time dynamical systems. In particular we shall show that strong 
anomalous diffusion may occur, i.e. 

((x(t) -x(0)) q ) (5) 

where v(q) is non-constant. In this case the probability distribution P(x,t) cannot be 
described in term of a unique scaling exponent, as in the Levy flight [14]] or in the random 
shear flow. 

Section 2 is devoted to a brief review of the known mechanisms for anomalous diffu- 
sion. In Sec. 3 we report detailed numerical results for the anomalous diffusion mechanism 
in 2d time-dependent incompressible velocity field, in the standard map (which can be 
considered a discrete time version of the Lagrangian motion in a 2d time periodic incom- 
pressible velocity field), and in an intermittent Id map. In Sec. 4 the reader can find some 
conjectures and discussions. 



2 A brief review on some aspects of the anomalous 
diffusion 

Anomalous diffusion occurs when some, or all, of the hypothesis of the central limit 
theorem break down. Practically, the system has to violate at least one of the two following 
conditions : 

a) finite variance of the velocity; 

b) fast enough decay of the auto-correlation function of Lagrangian velocities. 

This is, of course, an elementary remark. Nevertheless, in our opinion, a detailed dis- 
cussion of this point can be useful. In the present paper we use the term anomalous 
diffusion to indicate a non standard diffusion in the asymptotic limit of very long time. 
Sometimes in the literature the term anomalous is used also for long (but non asymptotic) 
transient behavior. This transient anomalous behavior can have practical relevance, e.g. 
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in geophysics. We do not discuss this point in details. Let us recall some well established 
results on the relation between diffusion phenomena and correlation function. It is easy 
to obtain the following formula : 

(( Xi (t) -x 4 (0)) 2 ) = f dh [ dt 2 {v i {x{t l ))v i {x{t 2 )))~2t [ drCair) (6) 
Jo Jo Jo 

where Cy-(r) is the correlation function of the Lagrangian velocity, namely v = x, 

C lJ (r) = (v t (x(r))v J (x(0))). (7) 

Eq. djD is the so-called Taylor relation ]15| for a test particle evolving according to the 
Langevin equation (|l|). 

For test particles whose dynamics is described by a map (discrete time), 

x(n + 1) = F(x(n)), (8) 

the Taylor formula (^j is still valid with the obvious changes : v(x(t)) i— > F(x(n)) — x(n) 
and J q Ch(t) i— > Cjj(0)/2 + $^j=i Ca(j). If both (v 2 ) < 00 and j dr Ch(t) < 00 then one 
has standard diffusion and the effective diffusion coefficients are 

1 f°° 
Di = lim —<(*,(*) - x,(0)) 2 ) = / dr C tl (r) (9) 

From the above remarks it stems that anomalous diffusion occurs only in two circum- 
stances: 

a) (v 2 ) = 00 

b) (v 2 ) < 00 and f* dr Cu{t) = 00, i.e. Cu(r) ~ r _/3 with (5 > 1, that means very 
strong correlations. 

2.1 Simple systems exhibiting anomalous diffusion 

The Levy flight (JTJJ belongs to the first case. Let us briefly discuss it for discrete one 
dimensional systems. The position x(t + 1) at the time t + 1 is obtained from x(t) as 
follows: 

x(t + l)=x(t)+U(t), (10) 

where U(t)'s are independent variables identically distributed according to a a-Levy- 
stable distribution, P a (U), with the following well-known properties : 

J dUe ikU P a (U) oc e- c|fc|Q and P a {U) ~ for | U |> 1 (11) 
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with < a < 2. An easy computation gives 

[ oo q > a. 

Note that (x 2 ) = oo for any a < 2, nevertheless one can consider the Levy flight as a sort 
of anomalous diffusion in the sense that x ~ t l ^ a ^> t 1 ^ 2 . 

In spite of the relevance of the a-Levy- stable distribution in probability theory, the 
Levy flight model, in our opinion, has a rather weak importance for physical systems. 
This is due to the very unrealistic property of infinite variance. 

An elegant way to overcome the above trouble is the introduction of a stochastic model, 
called Levy walk [|1(J, corresponding to ([!(]) but now U(t) is a random variable with 
nontrivial correlations. The velocity U(t) can assume the values ±uq and it maintains 
its value for a duration T which is a random variable with probability density i^{T). 
Practically, the particle moves with constant velocity for a certain time after which it 
changes direction and so on. In the Levy walk the origin of the possible anomaly is 
transferred to the correlation function of the Lagrangian velocity. Actually, taking 

ij){T) ~T~ (Q+1) , (13) 

one has standard diffusion if a > 2 while for a < 2 one has anomalous (super) diffusion : 

( 1/2 a > 2 

{x{t) 2 )~t 2v v=< (3-a)/2 Ka<2 (14) 
{ 1 a<\. 

Sometimes in the literature the Levy walk is (erroneously) called Levy flight. 



2.2 Nontrivial anomalous diffusion 

From the above considerations one realizes that it is not particularly difficult to build up 
ad hoc probabilistic models exhibiting anomalous diffusion. On the contrary it is much 
more difficult and much more interesting the understanding of the anomalous diffusion 
in nontrivial systems like the transport in incompressible velocity fields or deterministic 
maps, where the explicit probabilistic aspect (given for example by an external noise 
related to the molecular diffusivity) cannot play a relevant role. 



Avellaneda and Majda |T7|], (see also Avellaneda & Vergassola |TB| for the generaliza- 
tion to the time dependent case), obtained a very important and general result about the 
diffusion in an incompressible velocity field u(x). If the molecular diffusivity Do is non 
zero and the infrared contribution to the velocity field are weak enough, namely 

dk A «<»> I 2 ) < oo (1B) 



/ 
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where ( • ) indicates the time average and ii represents the Fourier transform of the 
velocity field, then one has standard diffusion, i.e. the diffusion coefficients D^s in (^) 
are finite. Therefore there exist two possible origins for the superdiffusion 

a) D > and, in order to violate the (|T5|), the velocity field with very long spatial 
correlation; 

b) Dq = and strong correlation between u(x(t)) and u(x(t + r)) at large r. 

For the sake of completeness we briefly remind one of the few nontrivial systems for which 
the presence of anomalous diffusion can be proven in a rigorous way. 
Consider a 2d random shear flow : 

u = (u(y) , 0) (f 6) 

where u(y) is a random function such that 

POD 

u(y) = dke iky u(k) (u(k)u(k f )) = S(k) 5(k-k f ), (17) 



S(k) is the spectrum and the average ( • ) is taken over the field realizations. Matheron 
& De Marsily || showed that the anomalous diffusion in the x-direction occurs if 

dk - = oo. 

On the contrary, if this integral is finite one has standard diffusion and 

Dfl=Da+ ±-i°° dk m. (18) 



Eq. (fL8|) is a well-known exact result already obtained by Zeldovich |T9[]. For the sake of 
simplicity we consider the case 

S{k) ~ W k^O. (19) 

If 7 > 1 then one has standard diffusion. On the contrary, if — 1 < 7 < 1 one has a 
superdiffusion, namely 

(|x(t)-*(0)| 2 >~^ ^ = ^>\ (2°) 

The condition f dk S(k) k~ 2 = 00 for the anomalous diffusion can be understood by 
means of the following simple physical argument. Note that J dk S(k) k~ 2 ~ (u 2 ) L 2 
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where L is the typical length of the function u(y), i.e. the typical distance between two 
sequent zeros of u(y). If (u 2 ) < oo and J dk S(k) k~ 2 < oo the diffusion process is 
basically similar to that one characterized by a velocity field given by a sequence of strips 
of size L and velocity ±a/ (u 2 ) and therefore the (|l8|) is nothing but the transversal Taylor 
diffusion p0[ in channels. The origin of the anomalous diffusion in this case is due to the 
fact that a test particle travels in a given direction for a very long time before it changes 
direction and so on. In some sense the Lagrangian motion in the random shear flow for 
— 1 < 7 < 1 is a nontrivial realization of a Levy walk. The generalization to the 3d case 
u = (u(y,z),0, 0) is straightforward. 

Another nontrivial system where the presence of anomalous diffusion has been rigor- 
ously proven is the Kraichnan's passive scalar model PTJ , where the velocity advecting the 



scalar field, is rapidly varying in time. For such a model, generally, the distance between 
pairs of particles tends to increase with the time elapsed but, occasionally, particles may 
come very close and stay so; this is the source of the anomalies in the scaling |22| |23| . 

If the inequality fll5|) holds then in order to have anomalous diffusion one needs Do = 
and a correlation function Ca(r) going to zero not too fast. These conditions are satisfied, 
e.g., in time periodic velocity fields whose Lagrangian phase space has a complicated self- 
similar structure of island and cantori J!]]. In this case super diffusion is essentially due to 
the almost trapping of the trajectories for arbitrarily long time, close to the cantori which 
are organized in complicated self-similar structures. 

In addition, there exist a strong numerical evidence (and sometimes also theoretical 
arguments) for superdiffusion in 1-d intermittent maps as well as in the standard map, 
which is a 2-d symplectic map and therefore it can be considered as the Poincare map of 
a Lagrangian motion in a time periodic velocity field. 

2.3 About the meaning of anomalous 

Up to now we used the term anomalous as synonymous of non-Gaussian. Let us consider 
now in more detail the cases with superdiffusion. There are two possibilities : 

a) weak anomalous diffusion when 

(| x(t) - x(0) \ q ) ~ t qv Vg > and v > - (21) 



b) strong anomalous diffusion when 



(| x(t)-x{0) \ q ) ~t qu ^ v{q)^cost u(2)>- (22) 
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where u(q) is a nondecreasing function of q. 
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In the weak anomalous diffusion one has for the probability distribution P(Ax, t) of 
Ax = x(t) — x(0) at time t, the following expression : 



P(Ax,t)=r u F(Axt- / ) (23) 

where in general the function F is different from the Gaussian one. The random shear 
flow is an example of the weak anomalous diffusion ||. 

Bouchaud at al. [24]] have shown that in the case of the random shear flow fll6|) -(|l8D 



with 7 = 

F{a) ~ e - c|a|4/3 for |a|>l. (24) 
There are very few numerical studies on the strong anomalous diffusion. Let us men- 



tion the work of Sneppen fa Jensen [25] on the motion of particles passively advected by 
the dynamical membranes (the authors use the term multidiffusion) and an 1-d intermit- 
tent map studied by Pikovsky |26| . 

Of course the strong anomalous diffusion is not compatible with the simple scaling 
fl2"3"|). Obviously in the presence of anomalous diffusion P(x,t) cannot obey to the usual 
Fickian diffusion Eq. (§). Therefore, a quite natural problem is to find the effective 
evolution equation for P(x, t) at large t. The answer to this question up to now as far as 
we know, is not well understood in the general case. 

For the Levy flight one can naively write 

d t P = -c (-8 2 ) a/2 P- (25) 

The above equation can be considered the proper one for the Levy flight in the sense 
that starting at t — with P(x, 0) = 5(x) at any t > the probability P(x,t) is given 
by th a-stable Levy function. Of course, if a = 2 one has the usual Fick equation and 
the Gaussian shape for P(x,t). For the more complicated anomalous diffusion there are 
many proposal in terms of fractional time, and/or spatial, derivatives. We do not enter 
in a detailed discussion of these works, the interested readers can see 



3 Numerical results 

In this section we report detailed numerical results for anomalous diffusion in 2-d time 
dependent velocity fields and maps. 
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3.1 Strong anomalous diffusion in a simple flow mimicking the 
Rayleigh— Benard convection 

We investigate here anomalous diffusion in a simple model mimicking the Rayleigh-Benard 
convection |28|, 2£|. Two-dimensional convection with rigid boundary conditions is de- 



scribed by the following stream function: 

i/;(x,y,t) = ^osin(x + Bsinut) siny , (26) 



where the periodicity of the cell is L = 2tt, and the even oscillatory instability |30 



is accounted for by the term B sin ut, representing the lateral oscillation of the rolls. 
The capability of the simple flow (|26|) to capture the essential features of the convection 
problem is discussed in Ref. f29|. 



It is worth noting that, at fixed B, the only relevant parameter controlling the diffusion 
process in the flow (p6|) is e = ujL 2 /ip Q , i.e., the ratio between the lateral roll oscillation 
frequency and the characteristic frequency of the passive scalar motion. Different regimes 
take place for different values of e. The two limiting cases e <C 1 and e ^> 1 have been 
investigated in Ref. J51J]. The regime e ~ 1 has recently been studied in Ref. [02J, in 



particular the synchronization between the lateral roll oscillation frequency (of the order 
of uj) and the characteristic frequencies of the scalar field motion (of the order of ipo/L 2 ) 
has been investigated. This mechanism leads to enhanced diffusion. 



We briefly recall the main results of Ref. [ 32 1 . It has been found that the synchroniza- 



tion between the circulation in the cells and their global oscillation is a very efficient way 



of jumping from cell to cell. This mechanism, similar to stochastic resonance |33|, [34|], 
makes the effective diffusivity as a function of the frequency to very structured. In the 
limit where the molecular diffusion vanishes, anomalous superdiffusion takes place for 
narrow windows of values of u around the peaks. 

In Fig. |I] the x-component, D^, of the eddy diffusivity versus the frequency ujL 2 /i/jq 
is shown for different values of Dq/i/jq. 

We anticipate that anomalous diffusion in the strong sense (according to the definition 
given in Sec. 2.2) takes place, making highly nontrivial the diffusion process only around 
the values of to corresponding to the peaks of enhanced diffusion. Specifically, let us 
focus our attention on the anomalous regime occurring for Dq = and uL 2 /ipo ~ 1.1. 
We have then integrated the Eq. (11) with Dq = and u obtained from the flow 
using a second-order Runge-Kutta scheme. In the following, averages are over different 
realizations and performed by uniformly distributing 10 6 particles in the basic periodic 
cell. The system evolution is computed up to times of the order of 10 4 L 2 /ip . 

In order to investigate the anomalous diffusion in the strong sense, a measure of 
([x(t) — x(0)] q ) for different g's has to be made. 
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Figure 1: The turbulent diffusivity D^/i/jq vs the frequency uL 2 /ipo for different values 
of the molecular diffusivity D /ip . D Q /ip Q = 3 x 1CT 3 (dotted curve); -Do/V'o = 1 x 
10~ 3 (broken curve); Dq/i/jq = 5 x 10~ 4 (full curve). 



To that end, we have performed a linear fit of ln([x(t) — x(0)] q ) versus lnt, the slope 
of which being the exponent qv(q) in (22). As an example, in Fig. the mean-squared 
displacement, as a function of the time, has been presented. Superdiffusive transport 
takes thus place for uj = 1.1 in the limit D — > 0, i.e.: 

([x(t) -x(O)f) oct 2i/(2) with i/(2) > 1/2 . (27) 

The presence of genuine anomalous diffusion is also supported by the behavior of the 
coefficient as a function of D . As suggested in Ref. J3S[, in the presence of genuine 



anomalous diffusion, the effective diffusivity must diverge and it is expected that Df ± ~ 
Dq 13 with > 0. The curve Dfi versus Dq is shown in Fig. |^. The data are well fitted 
by a straight line with slope (3 ~ 0.18, confirming the presence of anomalous diffusion at 

A) = o. 

Fits similar to that one presented in Fig. have been obtained also for other values 
of the order q, ranging between and 6. The results are summarized in Fig. f|, where the 
exponents qv(q) J s are shown as a function of q. 

Some remarks are in order. The curve, as a whole, displays a nonlinear behavior, 
the first clue of anomalous diffusion in the strong sense. Two linear regions are actu- 
ally present: the first one behaves up to q ~ 2, the second elsewhere. The two linear 
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Figure 2: The mean-squared displacement vs the time for the flow (p6|) with Dq = 0, and 
u = 1.1. Lengths and times are shown in units of L and L 2 /ipo, respectively. The best-fit 
(dashed) line corresponds to 2i/(2) = 1.3. 



regions are associated to two different mechanisms in the diffusion process, as suggested 
by the following simple considerations. For small g's, i.e. for the core of the probability 
distribution function F(Ax,t), only one exponent v\ = u(q) ~ 0.65 for q<2 fully char- 
acterizes the diffusion process even if anomalous. The typical, i.e. non rare, events obey 
a (weak) anomalous diffusion process, roughly speaking, one can say that at scale I the 
characteristic time r(Z) behaves like 

r(l)~l llv \ (28) 

On the other hand, the behavior qv(q) ~ q — 1.04 suggests that the large deviations are 
essentially associated to ballistic transport 

t{1) ~ I (29) 

basically due to the mechanism of synchronization between the circulation in the cells 
and their global oscillation. 

The character of strong anomalous diffusion, already evident from Fig. |4], is highlighted 
by the behavior of the probability density function (p.d.f.), P(Ax,t), of the displacement 
Ax = x(t) — x(0) at different times t's. As already pointed out in Sec. 2.3, weak and 
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Figure 3: The diffusion coefficient Df x as a function of Dq for the frequency of the roll 
oscillation uj = 1.1. The diffusivities are reported in units of ipQ. The best-fit (dashed) 
line has the slope —(3 = —0.18. 

strong anomalous diffusion can indeed carefully be discriminated by means of the sim- 
ple scaling law (23), which is not compatible with the strong anomalous diffusion. We 
have considered three different probability distribution functions for three different times, 
namely: Pi (Ax, ti), P2(Ax, £2) and P 3 (Ax, £3), with i 3 = 2t 2 = 4ii and ti, ti and £3 are 
well inside the scaling regions of the moments of the displacement. In Fig. [5] the three 
renormalized probability distribution functions P(Ax(t)/a), with a = exp(ln | Ax(t) |), 
at t — t\, and t 3 , respectively, are shown. 

Some remarks are in order. The first one concerns the tails of the three probability 
distribution functions which are clearly much higher than those of the Gaussian one. The 
strong character of the anomalous diffusion can be easily investigated by observing that 
if the scaling law (23) holds then the three functions rescaled probability distributions at 
different times have to be superimposed one on the other. As one can see from Fig. [5], this 
is not the case in our system. In fact, one satisfies (23) but only for small deviations, that 
means that only the cores of the probability distribution functions obey the above scaling 
law. This fact is not surprising being a direct consequence of the behavior shown for the 
scaling exponents gi/(g)'s. It is worth noting that the anomalous behavior disappears for 
a very small changes of the parameters, e.g. for u = 1.075 standard diffusion takes place. 
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Figure 4: The measured scaling exponents gz/(g)'s (joined by dot-dashed straight lines) 
of the moments of the displacement A function of the order q. The dashed line 

corresponds to 0.65 q while the dotted line corresponds to q — 1.04. 

3.2 Standard and anomalous diffusion in the standard map 

The standard map is a well-known symplectic map [[36] , |37j , given by: 



Jt+i = Jt + ^sin(27T^) , . 

0t+i = t + J t+ i mod 1. {6U> 

A short review of its phenomenology follows. For A = it is integrable 

Jt = Jo Vt /„, x 

e t = 6 + tJ modi; 1 J 

for all rational Jo one has periodic motion on resonant tori. As soon as A > the 
resonant tori disintegrate in pairs of elliptic and hyperbolic fixed points according to the 
Poincare-Birkhoff theorem [ 3l| . In the corresponding regions chaotic motions may occur. 



Still for small values of A the orbits remain trapped in invariant sets. This is due to 
the persistence under small perturbation of the nonresonant (KAM) tori where aperiodic 
motions take place and to the topological constraint imposed by the dimensionality of the 
system. In two dimensions the chaotic orbits generated by the disintegration of resonant 
tori are restricted to lie between two bounding KAM curves. Since the latter ones are 
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Figure 5: The normalized probability distribution function P(Ax(t)/cr) vs Ax = Ax/ a 
for the three times t\ = 500 (circles), ti = 2t\ (diamond) and £3 = It-i (squares). The 
dashed line represents the Gaussian function, shown for comparison. 

invariant and the area enclosed by them is conserved by symplecticity the chaotic regions 
are also invariant and no diffusion can take place |36| 

In order to observe diffusive behaviors one must consider values of A large enough so 
that all KAM tori which connect 9 = —0.5 to 9 = 0.5 disappear. The critical value of A 
has been numerically determined to be A c ~ 0.97 For A > A c the phase space can 
be pictorially described as a chaotic sea where islands of regular motions emerge. Such 
islands are related to the existence of accelerator modes i.e. stable periodic solutions of 
(|30"D specified by: 

Jt+p Jt = I f'?9"l 

AELi sin ( 2 ^) = 2ttZ 1 ' 

Each accelerator mode can be labelled with the integer pair (p,l). It is easy to check 
the stability for (1,1): the symplectic property implies that the Jacobian matrix of fl30D 
computed on a stable solution must give a pair of complex eigenvalues of magnitude one. 
The condition is satisfied for all A such that: 



2tt|/| < A < 2vr|/U/l + ( — ) (33) 



7T 
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Expression (p3|) implies that period one accelerator modes with different I cannot coexist. 
Furthermore the length of the interval of A values such that (1, 1) is stable decreases for 
large / as: 

AA ~ jjj (34) 

If only one stable island exists no sensitive contribution to the diffusion properties of the 
chaotic orbits is expected. A trajectory that arrives in the neighbourhood of a stable island 
sticks to it for a finite amount of time before getting away. On average the contribution 



to the encircling orbit is zero. Typically one observes normal diffusion |39| : 

([J t - J ] 2 ) ~2Dt t>l (35) 

where the average is over the initial conditions, the diffusion coefficient D can be estimated 
with a crude random phase approximation as: 

D,D QL = - J t ) 2 ) = i (A) 2 1 „ *? {Mf) = AL. (36) 

In Fig. |6] the standard diffusion behavior of the exponent q v(q) as a function of q is 
shown for A = 3.86 i.e. in the absence of any stable islands whereas. The same standard 



diffusion occurs for A = 10.492927 when one stable island exists []37[,||- For both the 
cases standard diffusion occurs. 

The coexistence of many accelerator modes has relevant consequences on the diffusive 
behavior of the system. For values of A corresponding to such a situation, the variance 
goes asymptotically as 

([Jt ~ Jo?) ~ t 2v{2) t > 1 with i/(2) ^ 1/2. (37) 

The sticking of the chaotic orbits to stable islands leads to the appearance of blocks of 
long range correlation in the sequences of the Jt variable. This behavior is qualitatively 
similar to the features of the Levy walk. 

For the case corresponding to A = 6.9115, where the fundamental accelerator mode 
(1, 1) coexists with the two orbits with period three (3, ±3), we have found 2z/(2) ~ 1.32 
in agreement with previous works ||. 

On the other hand, in order to investigate the presence of strong anomalous diffusion, 
the knowledge of z/(2) does not characterize completely the diffusion properties of the 
orbits. From the numerical results (see Fig. ^), one can extrapolate with fair agreement 

/ x f 0.59g q < 2 , ooN 
q ^ = \q-0m q>2. ^ 
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Figure 6: The diffusion exponent qv{q) as a function of q for A = 3.86. The results have 
been obtained by averaging over 10 6 trajectories with J = 0. 
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Figure 7: The diffusion exponent qv{q) as a function of q for A = 6.9115. The averages 
are over a sample of 10 6 trajectories. The dashed line corresponds to 0.59 q while the 
dotted line corresponds to q — 0.69. 
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Figure 8: The rescaled probability distribution function P(AJ(t)/a) after N = 4 x 10 6 
iterations, for A = 6.9115 and t\ = 10 2 (circles), £2 = 10 3 (squares), £4 = 10 4 (diamonds) 
with a(t) = exp{< In |Aj t | >}. 

The second moment u(2) ~ 0.65 falls in a crossover region between the two behaviors. 
The feature ( |38D suggests the presence of two mechanisms for the typical events and the 
rare (ballistic) ones as already discussed for the 2-d time dependent flow in the previous 
subsection. 

Such a conjecture can be tested by a direct inspection of the probability distribution 
function. In Fig. |] the rescaled probability distribution function of AJ t = 5J t /cx with 
a = exp{< In |AJ t | >} for different times t is shown. 

A further natural question to be addressed is how generic is the strong anomalous 
behavior observed for A = 6.9115 in the standard map. To inquire this point we have 
studied the diffusion properties for (a), A = 6.8115, (6) A = 7.0115 and (c) A = 6.4717. 

The rationale of (a) and (b) is to check the stability of the result (]38f) with respect to 
small changes of A. In the case (a) strong anomalous diffusion is observed. The behavior 
of the moments practically coincides with that one observed for A = 6.9115. For (b) A = 
7.0115 we observe normal diffusion: the fact has to be related to the destabilisation of 
the pair (3, ±3). Finally (c) corresponds to the coexistence with (1, 1) of two other stable 
islands now generated by the pair of period five accelerator modes (5, ±5). Also in this 
case there is fair evidence of strong anomalous diffusion although now the crossover region 
between the two regimes is more extended. 
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At the end of this subsection we briefly discuss some results recently obtained in 
4"2L E3[ where the transport in 2 — d symplectic maps (practically a modification of the 



standard map) is considered. The authors find that the anomalous diffusion is rather 
rare. This is particularly clear in Figs. 5 and 6 of ||2| where the exponent z/(2) for the 
kicked Harper map is shown to be always equal to 1/2 apart very narrow regions in the 
control parameter space. 

3.3 Anomalous diffusion in 1-d intermittent maps 

Strong anomalous diffusion may also occur in simpler dynamical systems of the form 

x t+ i = f(x t ) (3g) 

yt+i = yt + xt 

where f(x) is a symmetrical Lorenz type map on the interval [—1, 1]. Diffusion may then 
appear in the variable y t . As an example we reconsidered a map previously studied by 
Pikovsky |]26[| . The function f(x) is implicitly specified, Vz > 1, for x g]0, 1] by 

x -{ £[1 + /(*)]* f °r 0<*<^ 

\ f(x) + £[1 -/(*)]* for i<x<l [W) 

with the further requirement of monotonicity. The extension to the interval [— 1,0[ is 
obtained by imposing 

f(-x) = -f(x) (41) 

The resulting shape of f(x) is plotted in Fig. ||. 

Thank to the implicit definition fl4~0|), by a straightforward computation, it is easy to 
verify that the map has a uniform density invariant distribution. 

As discussed in Sec. [| anomalous diffusion follows from power-law behavior of the 
correlation function of the increments of the diffusing variable. Here the latter property 
stems from the intermittency of the dynamics. The mapping (|4*0| ) induces in the interval 
[— 1, 1] two qualitatively different regions. It is not restrictive to focus on [0, 1] 

In the neighbourhood of x = 1 the orbits perform laminar motions. The basic features 
of the dynamics are captured || by the polynomial behavior around the unstable fixed 
point. 

f{x) = x + liz^ll + o(( i _ x y) (42) 

^From (|4*2] ) it follows that the typical duration of the laminar phase for an orbit starting 
from x ~ 1 is T x ~ (1 — x)^^^. Since the invariant measure is uniform a reasonable 
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Figure 9: f(x) vs x from Eq. (j42[) for z = 2.5 (full line) and 1.5 (dashed line). 



inference j40| is that the probability distribution of having laminar motion of duration t 
scales as 

Plam{t) ~ t-^. (43) 

A qualitatively different picture comes from the set [0, j-]. There the mapping has 
the explicit representation: 

f{x) = [(2z\x\)t-l]. (44) 

In this second region the motion is chaotic and the increments Ax t = x t +i —x t decorrelate 
after few steps. The coexistence of the two regimes prevents the dynamics from having a 
typical timescale so that the increments correlation function obeys to a power law 

C(r) =< Ax t+T Ax t >~ P lam (t) ~ r 1 '^ (45) 

as it has been shown in Anomalous diffusion is therefore expected for z > 2. 

We studied the diffusion behavior of the dynamics for four values of z. 

For z = 1.5 the integral of the correlation fl4"5]) is convergent therefore the second 
moment grows linearly with time. We observe standard scaling for all the moments of 
power q < 3 as shown in Fig. |10 . 
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Figure 10: The measured scaling exponents qu(qys of the moments of the displacement 
Axt, as a function of the order q for z = 1.5. The dashed line 0.5 q while the dotted line 
corresponds to q — 1.50. 

Let us note that in this case, since 2/(2) = 1/2, the diffusion seems to be standard. 
Actually, because of the nonlinear behavior of qv(q), one can say that this is an example 
of strong anomalous diffusion. 

For z = 2 the mapping coincides with flH| ) in all the interval [0,1]. This allows the 
use of a faster iteration algorithm. The value z = 2 is also critical for the appearance 
of anomalous diffusion, in the sense u(2) > 1/2, due to the divergence of X^li C(r). 
Actually the plot of the diffusion exponent qv(q) versus q gives a fair evidence for strong 
anomalous diffusion. In Fig. [FT] the moments larger than the fourth are in agreement 
with a ballistic dynamics qv(q) ~ q — 0.54 while for moments with q < 2 it results 
q u(q) ~ 0.61 q. 

The non Gaussian behavior of the core of the probability distribution function appears 
also from the direct inspection of the distribution as shown in Fig. [12] where the probability 
distributions functions calculated for different times and rescaled by a are compared with 
the normal distribution. 



In Fig. [13] and 14] q u(q) versus q is reported for z = 2.5 and z = 3.0. 



As in [E9] we can infer from the previous analysis the following form of the PDF of 
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Figure 11: The same as in Fig. [T(J for z = 2.0. The dashed line corresponds to 0.61g while 
the dotted line corresponds to q — 1.03. 
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Figure 12: The rescaled probability density as in Fig. || but for the Pikovsky map when 
z = 2 in Eq. (M). 
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Figure 13: The same as Fig. |10| but for z = 2.5. The dashed line corresponds to 0.73g 
while the dotted line corresponds to q — 0.54. 
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Figure 14: The same as Fig. |10| but for z = 3.0. The dashed line corresponds to 0.78g 
while the dotted line corresponds to q — 0.44. 
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the diffusing variable y t for z > 2: 



P(y, t) = [l- P lam (t)]r Ucore F core (Jr»—) + P lam (t)[5(y - t) + 5{y + t)}. 

(46) 

A consequence of (fPf) is 

g%) = ( ^ ore n , g< J (47) 
^ w/ [ g - l/(z -1) q> 2 v 7 

The predictions of (^) are in qualitative agreement with the numerical experiments. For 
z < 2 the core of the probability distribution functions becomes Gaussian. 



4 Conjectures and remarks 

We have given rather clear numerical evidence of the presence of strong anomalous dif- 
fusion generated by nontrivial dynamics. Let us do some remarks and conjectures. As 
first we note (Sec. 3) that the anomalous diffusion is, in some sense, rare. For example 
in the flow mimicking the Rayleigh-Benard convection (Sec. 3.1) the anomalous diffusion 
occurs only for some precise values of the parameter u. Something similar, even if in 
a weaker way, happens for the standard map (Sec. 3.2) or other 2-d symplectic maps 



42], [4J|. Even though for the 1-d intermittent maps anomalous diffusion could seem a 
generic feature this is not the case. Actually the origin of the super diffusion is in the 
very peculiar property (marginal instability) of the fixed point (see Eq. (|4~2|) ) destroyed 
by any generic modifications of the map. Asymptotic anomalous diffusion, at least that 
one generated by non ad hoc dynamical systems, seems to be a nongeneric property i.e. 
it disappears as soon as a small perturbation 0(e) in the evolution law is introduced. 
However transient anomalous diffusion appears in the perturbed system up to a crossover 
time t* ~ e~ a while for larger time standard diffusion takes place with a rather large 
diffusion coefficient. Therefore one can say that, even if asymptotic anomalous diffusion 



is very rare, its ghost is visible |35 



In all the cases here considered strong anomalous diffusion has been found to corre- 
spond to a rather simple shape of the z/(g): 



av( q )^\ Uiq q<<lc (48) 
* yH ' \ q - c q > q c . v ' 

The above behavior suggests the existence of two different mechanisms in the diffusion 
process. The typical (i.e. non rare) events obey a (weak) anomalous diffusion process 
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characterized by the exponent u%, i.e. the characteristic time r(Z) at the scale / behaves 
like r(l) ~ l 1 /" 1 , but the large deviations are essentially associated to a ballistic transport 
r(Z) ~ /. In the flow mimicking the Rayleigh-Benard convection the ballistic mechanism 
is due to the synchronization between the circulation in the cells and their oscillations. In 
the standard map it arises from the accelerator modes while in the 1-d intermittent maps 
is due to the marginally unstable fixed point. 

We do not claim of course that the above shape fl4lf ) is universal. It is rather easy, in 
fact, to build up a toy system for which u(q) can be any function with the only constraint 
that q u(q) is a convex function. Following Elliott et al. HI] we consider the velocity field 



u=(u(y),w) (49) 

where w is a constant and u(y) is a given quenched random field. Neglecting the molecular 
diffusion the Lagrangian evolution is governed by the following equations : 

* = <y) m 
y = w. 

Now, taking for the sake of simplicity w = 1, and y(0) = the following equation holds 

x(t) - x(0) = [ dt' u{t') (51) 



Jo 

which trivially relates the Lagrangian and the Eulerian statistical properties of the system. 
Thus, from Eq. ( |5lD it is not difficult understand that adopting a suitable multiaffine 
process [^, [RJ for u(y) one can find the relation defining strong anomalous diffusion 

(| x(t) - x(0) H ~ t qu{q) v{q) + cost i/(2) > - (52) 

with any proper shape for v{q). A fast and efficient algorithm to generate a multiaffine 
process is done in |47| . 



A very interesting open problem is the determination of the effective diffusion equation 
both at long times and large scales for a system with strong anomalous diffusion. We stress 
the fact that if u(q) ^ const the probability distribution P(Ax, t) cannot obey a simple 
generalization of the standard diffusion equation. For sure a P(Ax, t) corresponding to a 
v(q) such that ( |52"D holds cannot be obtained from generalized linear diffusion equation 
involving fractional spatial derivative as Eq. (|25|). Even other recently proposed linear 
diffusion equations with fractional derivatives both in space and in time are not able to 
produce such nontrivial behavior |27 |. 
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